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1. Main text 



In the last two decades, several hundred planets have been detected beyond our Solar-system. 
Most of these extra-solar planets orbit sun-like stars. A small number have been detected around 
stars that are in their late evolutionary state, such as Red Giant Branch (RGB) stars and pulsars. 
The phase directly after the RGB stage, the Horizontal Branch (HB), however, is still unexplored; 
therefore, there is no empirical evidence for whether close-in planets, i.e., those with semi-major 
axes less than 0.1 AU, survive the giant phase of their host stars. 

Besides its evolutionary stage, a star's chemical composi tion appears to be a major indicator 



Gonzales 



(1997), showed that 



of its probability for hosting a planet. Previous studies, e.g., 
main-sequence (MS) stars that host giant planets are metal-rich. This finding is supported by the 
large exoplanet search surveys around MS stars reporting a c onnection between planet frequency 



and metallicity (lSantos et al 



2004; 



Valenti & Fischer 



main- sequence stars finding no evidence for Jovian planets (|Sozzetti et al. 



2005|), and a survey of 1 60 metal-poor 



2009) 



Until now, only very few planets have been detected around stars with metallicities as low 
as [Fe/H]= — 1, i.e. 10% of the sun's metallicity. The detection of PSR B1620 b, a Jovian planet 



orbiting a pulsar in the core of the metal-poor globular c 



however, that planets may form around metal-poor stars (IFord et al. 



uster M4 ( 



Fe/H 



2000; 



=-1.2), suggests, 



Sigurdsson et al. 



2003(), although the f ormation mechan ism of this particular planet might be linked to the dense 



cluster environment (Bee r et al. 



2004). 



We used the Fibre-fed Extended Range Optical Spectrograph (FEROS), a high-resolution 
spectrograph (R = 48, 000) attached to the 2.2 meter Max-Planck Gesellschaft/European Southern 
Observatory (MPG/ESO) telescopej, to observe the star HIP 13044. This star is classified as 



'The observations of HIP 13044 were carried out from September 2009 until July 2010. The 
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a red HB (RHB) star 



Chiba & Beers 



2000; 



' Fig. 1) and its meta 



Carney et al. 



2008b; 



content is TFe/H J 



mean 



-2.09 (lBeersetal.ll 19901 : 



Roederer et al. Il2010h . i.e. about 1% that of the Sun. 



So far, HIP 13044 is not known as a binary system. Detailed stellar parameters can be found in 
Supporting Online Material, Section I2TT1 

Previous radial velocity (RV) measurements of HIP 13044 showed a systematic velocity 



of about 300 km s with respect to the Sun, indicating that the star belongs to t he stellar ha 



(Car ney et al. 



1986j) . Indeed, the star has been connected to the Helmi stream (|Helmi et al 



o 



1999), a group of stars that share similar orbital parameters that stand apart from those of the 
bulk of other stars in the solar neighborhood. The Helmi stream stars move on prograde eccentric 
orbits (-Rperi ~ 7 kpc, R apo ~ 16 kpc) that reach distances up to |z| max ~ 13 kpc above and below 
the galactic plane. From t hat, it has been concluded that these stars were once bound to a satellite 



galaxy o f the Milky Way dHelmi et al. 



Ga ago dKepley et al. 



2007) 



19991 : 



Chiba & Beers 



2000) that was tidally disrupted 6-9 



The variation of the RV between our observations at different epochs has a semi- amplitude of 
120 m s -] c| (Fig. 2). The Generalized Lomb Scargle (GLS) periodogram Jzechmeister & Kiirster 



2009) reveals a significant RV periodicity at P = 16.2 days with a False Ala rm Probability of 



5.5 x 10 6 . Additional analysis, using a Bayesian algorithm ( 



Gregory 



2005), yields a similar 



period around 16 days. Such RV variation can be induced by an unseen orbiting companion, 



spectrograph covers a wavelength range from 350 to 920 nm (|Kaufer et al 



2000). To measure 



the RV values of HIP 13044 we used a cross-correlation technique, where the stellar spectrum 
is cross-correlated with a numerical template (mask) designed for stars of the spectral type F0 
(Section |Z2]>. 

2 In order to search for periodic variations, we used periodogram analysis techniques, which are 
capable of treating missing values and unevenly spaced time points. 
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by moving/rotating surface inhomogeneities or by non-radial stellar pulsations. Exploring both 
stellar rotational modulation and pulsations is critical when probing the presence of a planetary 
companion, because they can produce a similar or even the same RV variation, mimicking a 
Keplerian motion. 



A well-established technique t o detect ste 



line profile asymmetry or bisector (|Gray et al 



ar ro tational modulation is to investigate the 



2008b and Ca II lines (Section Q. Surface 



inhomogeneities, such as starspots and large granulation cells, produce asymmetry in the spectral 
line profiles. When a spotted star rotates, the barycenter of the line profiles moves periodically and 
the variation can mimic a RV variation caused by an orbiting companion. Instead of measuring 
the bisectors , one can equi valently use the bisector velocity spans (BVS) to search for rotational 



modulation (Hatzes 



19961) . Adopting this technique, we have measured BVS from the stellar 
spectra. There is only a weak correlation between BVS and RV (correlation coefficient =—0.13), 
but the BVS variation shows a clear periodicity with P =5.02 days (Section [2.3 .11) . No period 
around 16 days is found in the BVS variation. 

In addition to the BVS analysis, we investigated the variation of the Ca II A849.8 line, 
which is one of the Ca II infrared triplet lines. From the observed Ca II 849.8 equivalent- width 
variations we computed a mean period of 6.05 days (Section [2.3.11) . which is in the same order of 
the period of the BVS variation. We adopted the mean period of both stellar activity indicators, 
-Prot = 5.53 ± 0.73 days, as the stellar rotation period of HIP 13044 and then calculated the 
inclination angle of the stella r rotation axis, which f ollows from P ro t/ sini = 2-KR*/vsvai. With 



a stellar radius iL 



6.7 R ICarnev et al. 



(2008a) and our adopted value for the projected 
rotational velocity, v sinz=10.25 kms -1 , which was derived from the observed line broadening 
(Section [27Tb . we obtained an inclination angle i = 9.7 ± 1.3 deg. Thus, the real stellar rotation 
velocity is ~62 kms -1 , which is typical for an early F type MS-star but relatively high for HB 
stars. 
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An explanation for this high rotation velocity is the assumption that HIP 13044 has engulfed 



its close-in p 



(Soker 



1998 



anets during the red 



Levrard et al. 



2009 



giant phase. Infalling planets are able to spin-up their host star 



Carlberg et al. 



2009), and this mech anism has been sugg ested 



to explain the high v sin i values observed for many RGB and HB stars (|Carney et al. 



2003). 



We observed vari ations of HIP 1304 4 in the photometric data from the Hipparcos satellite 



(|Perryman et al. 



19971) and SuperWASP (IPollacco et al 



20061) (Section [23^1) . While the 



Hipparcos data shows only a marginal significant periodicity of 7.1 hours (FAP=1.8%), the 
SuperWASP data shows few intra-day periodicities with FAP~1% and two significant periodicities 
at 1.39 (FAP=5 x 1CT 4 ) and 3.53 days (FAP=2 x 10~ 4 ). These two periods, however, are most 
likely harmonic to each other (1.4 -1 + 3.5 -1 = 1 ). It is expected tha t HIP 13044 oscillates only 



at pulsationally unstable overtones of high order (Xiong et al. 



star in the metal-poor globular cluster NGC 6397 (IStello & Gilliland 



predictions (IXiong et al 



1998). Observations of one RHB 



2009) as well as theoretical 



1998|) set these periods in the range of a few hours to a day or so. No 
clear theoretical predictions for a star with parameters similar to HIP 13044 exist, hence it is 
possible that some high-order oscillations are able to explain the 1.4 or 3.5 day signal. What is 
important, however, is that there is no signal of a period around 16.2 d in the photometric data. 

The arguments above show that neither stellar rotational modulation nor pulsations are 
plausible sources of the observed periodic RV variation. Therefore, the best explanation for the 
~16 days period is the presence of an unseen companion. We computed its orbital solution 
(Table 1). Its minimum mass lies securely in the planetary mass domain, even with a plausible 
sin % uncertainty. With an eccentricity of 0.25 and a semi-major axis of 0.116 AU, the planet 
has a periastron distance of about 0.087 AU which is ^2.8 times the present stellar radius. The 
periastron is ~0.06 AU away from the stellar surface. 



Because a large number of known exoplanets have orbital semi-major axes between 0.01 and 
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Table 1: Orbital paramete rs of HIP 13044 b 

P 16.2 ± 0.3 days 

K x 119.9 ±9.8 ms" 

e 0.25 ± 0.05 





219.8 ± 1.8 


deg 


JD - 2450000 


5109.78 ±0.02 


days 


x 2 


32.35 


ms~ ] 


rms 


50.86 


ms~ ] 


rri2sini 


1.25 ±0.05 


M Jup 


a 


0.116 ±0.01 


AU 



0.06 AU, the distance between the periastron and the star HIP 13044 itself is not unusual. The 
non-circular orbit (e =0.25), however, is not expected for a close-in giant planet around a post 
RGB star. 



In the case of HIP 13044, the original orbit could have been disturbed o r changed du ring 



the evolution of the star-planet-system, in particular during the RGB phase (|Soker 



1998). 



Interestingly, the orbital period of HIP 13044 b is close to three times the stellar rotation period. 
There are a number of known planetary systems which also have such a "coupling" between the 
stellar rotation and orbital periods, e.g. Tau Boo (1:1), HD 168433 (1:2), HD 90156 (1:2) and HD 
93083 (1:3). Such p lanetary systems are particularly interesting to study star-planet interactions 



(Shkolnik et al. 



2008) 



So far, there are only very few planet or brown dwarf detecti ons aroun d post RGB stars 



besides the pulsa r planets, namely V391 Peg (ISilvotti et al. 1120071) . HW Vir (ILee et al. 1120091) 



and HD 149382 (Geieret al 



2009) (Fig. 1). These are, however, substellar companions around 
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subdwarf-B or Extreme Horizontal Branch (EHB) stars, i.e., the nature of their host stars 
differs from that of HIP 13044, an RHB star. Contrary to RGB star s , such as G and K gia nts 



(Wal ker et al. 



subgiants (e.g. 



1992 



Hatzes & Cochran 



Johnson et al. 



1993 



Setiawan et al. 



2004; 



Dollinger et al. 



2007) and 



2010), HB stars have not been yet extensively surveyed for planets. 



While at least 150 main-sequence stars are known to bear close-in (a = 0.1 AU) giant 
planets, so far no such planets have been reported around RGB stars. A possible explanation is 
that the inner planets have been engulfed by the star when the stellar atmosphere expanded during 



the giant phase. The s urvival of HIP 1304^ 



certain circumstances (Livio & Soker 



1983 



~> during that phase is theoretical 



Soker 



1998 



Bear & Soker 



y possible under 



2010() . It is also possible 



that the planet's orbit decayed through tidal interaction with the stellar envelope. However, a 
prerequisite to survival is then that the mass loss of the star stops before the planet would have 



been evaporated or accreted. Assuming asymmetric mass loss, veloci ty kic 



cs could have increased 



2007). The same could 



the eccentricity of HIP 13044 b to its current, somewhat high value (IHeyl 
be achieved by interaction with a third body in the system. 

Interestingly, a survey to charact erize the multiplicity of EHB stars showed that more than 



60% of the sample are close binaries (M axted et al 



2001). Their orbital radii are much smaller 



than the stellar radius in the RGB phase. This could be explained by the high friction in the 
interstellar medium, which would move a distant companion towards the primary. Such spiral-in 
mechanism could also take place in the RGB-to-RHB transition phase. Similar to the binary 
case, a distant giant planet in the RGB phase can move towards the primary into a smaller orbit. 
Consequently, the resulting close orbiting planets will be engulfed when the stellar envelope 
expands again in the next giant phase. HIP 13044 b could be a planet that is just about to be 
engulfed by its star. 



HIP 13044, with a mean metallicity estimate of [Fe/H]=-2.1, is far more metal-poor than 
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any previously known exoplanet hosting star (Fig. 3). For the existing theories of hot giant 
planet formation, metallicity is an import ant parameter: in particular, it is fundamental for the 



core- accretion planet formation model dlda & Lin 



20051) . It might be that initially, in the planet 



formation phase, HIP 13044 had a higher metallicity, and that during its subsequent evolution, it 
lost its heavier elements. For example, during the giant phase, heavy elements could have had bee n 



incorp orated into dust grains and then separated from the star's atmosphere (|Mathis & Lamers 



1992|) . However, given the star's membership to the Helmi stream, in which the most metal-rich 



2009J), we do not expect its initial 



sub-dwarfs known so far have [Fe/H]~ —1.5 (|Klement et al. 
Fe abundance to exceed this value. 

Finally, as a member of the Helmi stream, HIP 13044 most probably has an extragalactic 
origin. This implies that its history is likely different from those of the majority of known 
planet-hosting stars. HIP 13044 was probably attracted to the Milky Way several Ga ago. Before 
that, it could have had belonged to a satellite galaxy o f the Milky Way similar to Fornax or the 



Sagittarius dwarf spheroidal galaxy (|Helmi et al 



1999). Because of the long galactic relaxation 



timescale, it is extremely unlikely that HIP 13044 b joined its host star through exchange with 
some Milky Way star, after the former had been tidally stripped. The planet HIP 13044 b could 
thus have a non- Galactic origin. 
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Fig. 1. — Location of HIP 13044 in a My vs. B — V color-magnitude diagram (CMD) shown as 
a black dot superimposed to the CMD of Messier 3 (grey open circles) based on the photometry 



by 



Buonanno et al. 



(I1994|) . Apparent magnitudes have be en converted to absolute magnitudes by 



Harris 



(I1996|) . The gap separating the 



considering the distance modulus and extinction given by 
blue and red Horizontal Branch (HB) is due to RR Lyrae instability strip. The CMD location of 
HIP 13044 implies that it is a core He-burning star, located a the blue edge of the RHB. Further 



candidates for post RGB stars 



( Silvottietal. 



2007 



hosting planets/brown dw arf, V391 Peg, HW Vir and HD 149382 



Lee et al. 



2009 



Geier et al. 



2009) are displayed as crosses. 



- 10- 



4.4 



4.2 



O 
O 

CO 



4 - 



HIP 13044 



1 



5100 



5200 5300 
Time (JD-2450000) 







0.5 



1.5 



Phase 



(a) 




5400 




Fig. 2. — (Upper panel) RV variation of HIP 13044. The RV values have been computed from 
the mean RVs of 20 usable echelle orders of the individual spectrum. The error bars have been 
calculated from the standard error of the mean RV of each order. (Lower panel) RV variation 
phase-folded with P =16.2 days. 
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Metallicity [Fe/H] 

Fig. 3. — Distribution of the metallicity [Fe/H] of planet-hosting stars. 
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2. Supporting online material 

This supplementary online material is part of the article with the title "A Giant Planet Around 
a Metal-poor Star of Extragalactic Origin". Here we give additional information about the stellar 
parameters of HIP 13044 (Section |27TT > and more detailed descriptions of the observations as well as the 
stellar radial velocity (RV) computations (Section I2721 ). As reported in the manuscript, the RV variation 
has a period of 16.2 days and semi-amplitude of 120 ms _1 which most likely originate from an unseen 
planetary companion. In the following sections we describe briefly the technique that we used to investigate 
the stellar activity (Section |23T >. We analyzed the line profile asymmetry (bisector) and Ca II A849.8 line, 
which are excellent stellar activity indicators. Additionaly, we also present our analysis of the Hipparcos 
and SuperWASP photometric data. 



2.1. Stellar parameters 



HIP 13044 (CD-36 1052, CPD-36 287) is a star of spectral type F2 (SIMBAD). It ha s a visual 



magnitude of my =9.94 mag (Hipparcos) and distance of ~700 pc ([Carney et al 



1986 



2008a) 



HIP 13044 resides in the red part of the Horizontal Branch (RHB), which is separated from the 
blue and extreme Horizontal Branch (BHB and EHB) stars by the RR Lyrae instability strip. 



The Hipparcos archive data (IPerryman et al 



19971) revealed variations in the photometry 



and astrometry. From the Hipparcos intermediate astrometry data, we calculated the rms of 
the difference between the observed great-circle abscissa and the solution computed with the 
reference astrometric parameters to be 9 x 10~ 3 arcsec. While the Hipparcos photometry data 
shows only a marginal significant periodicity at P = 7.1 hours, photometric observations with 



SuperWASP (IPollacco et al. 



2006) reveal some periodicities with timescales of several hours up 



to few days (Section r2.3.2l) . These signals might correspond to the stellar pulsations. 
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As a member of the Helmi stellar stream (IHelmi et al. Ill999h . HIP 13044 shares the property 
of other stream members, like the low iron abundances ([Fe/H]= —1.8 for 33 stream members) 



(Kep lev et al. 



(Roederer et al 



Roederer et al. 



2007; Klement et al. 



2009) and a chemical similarity to typical inner halo stars 



2010). An extensive abundance analysis of HIP 13044 has been also presented in 



(20K 



Table. SI 
Stellar parameters of HIP 13044 



Parameter 


Value 


Unit 


Reference 


Spectral type 


F2 




SIMBAD 


my 


9.94 


mag 


Hipparcos 


distance 


701 ±20 


pc 


Beers etal. ( 1990): Carnev et al. (1986. 2008a) 


T e fi 


6025 ±63 


K 


Carnevetal. (2008a): Roederer et al. (2010) 


i?* 


6.7±0.3 


R 


Carnev etal. (2008a) 


log g 


2.69±0.3 




Carnevetal. (2008a) 


m 


0.8±0.1 


M 


this work 


[Fe/H] 


-2.09±0.26 


[Fe/H] 


Chiba & Beers (2000): Carnev et al. (2008b): Roederer et al. 


v sini 


8.8±0.8 


kms -1 


Carnevetal. (2008a) 




11.7±1.0 


kms -1 


this work 



Table SI summarizes the fundamental stellar parameters, such as the effective temperature 
T eff , surface gravity log ;g, stellar radius R*, stellar mass and metallicity [Fe/H], from the 
literature and our spectroscopic measurements. The stellar mass has been inferred from the 
knowledge of th e stellar rad i us R* a nd surface gravity log g. We calculated a stellar mass of 



0.8±0.1M„. 



Carney et al. I (I2008ah have measured vsini = 8.8 kms , whereas we obtained 



v sin i = 11.7 k m s 1 . The discrepanc y between the two results is probably caused by the different 



methods used. ICarney et al. 



(2008a) used a Fourier transform method, whereas we used a 
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cross-correlation technique to measure the v sini (|Setiawan et al 
vsini = 10.25 ±2.1 kms -1 which is the mean value of both measurements 



2004). For this work we adopt 



2000). The first fiber 



2.2. Radial velocity measurements 

The observations of HIP 13044 have been performed with FEROS, a high-resolution 
spectrograph (R = 48 000) attached to the 2.2 meter MPG/ESO telescope, located at the ESO La 
Silla observatory. The spectrograph is equiped with two fi bers and has a wav elength coverage 
from 350 to 920 nm that is divided into 39 echelle orders (|Kaufer et al. 
is used to record the stellar spectrum, whereas the second fiber can be used for a simultaneous 
calibration with a ThArNe lamp, in which the instrumental velocity drift during the night can be 
corrected. 

The raw spectroscopic data have been reduced with an online data reduction pipeline that 
produced 39 one-dimensional spectra of each fiber. The stellar RVs have been computed by using 
a cross-correlation technique, where the stellar spectrum is cross-correlated with a numerical 
template (mask ) . The data reduction and RV computation methods have been described in 



Setiawan et al. 



(|2003r) . From our long-term RV surveys (2003 to 2010) with FEROS, an accuracy 
in RV of 11 ms" 1 has been measured for a spectroscopic RV standard star HD 10700 (spectral 
type K8V). 

However, this accuracy cannot be achieved for stars with fewer absorption lines due to higher 
effective temperature and for fast rotating stars, since the spectral lines are broader. We obtained 
a typical accuracy of ~ 50 m s _1 for the individual RV measurements of HIP 13044. In total, 37 
high-resolution stellar spectra have been obtained with FEROS. However, one spectrum could 
not be used for the RV analysis due to a calibration problem. Thus, we use 36 spectra for our 
spectroscopic study. For HIP 13044 (spectral type F2) we measured the stellar RVs by computing 
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the cross-correlations of the stellar spectra with a mask, designed for stars of the spectral type 
FO. In addition, we also created a special cross-correlation mask from the mean spectrum of 
HIP 13044 spectra. This mask, called "HIP 13044 temp", contains 550 selected spectral lines. To 
find out if there is inconsistency in the characteristics, such as period and amplitude, of the RV 
variation, as a result of using two different masks, we cross-correlated the spectra of HIP 13044 
with this special mask HIP13044-temp. Except the zero-point offsets, we found no discrepancies 
in the period and amplitude of the RV variation when using either the FO or the HIP13044-temp 
masks. Thus, in this work we present only the RV measurements with the standard FO mask. 

The RV values of HIP 13044 have been computed from the mean RVs of 20 usable echelle 
orders. Table S2 gives the values of the RV mea surements. A periodogram analysis by using a 



Lomb-Scargle (LS) periodogram ( 



Scargle 



19821) shows a significant peak at P ~16 days with 



a False Alarm Probability (FAP LS ) of 2 x 10 3 . In addition, we also applied the Generalized 
Lomb-Scargle (GLS) periodogram (|Zechmeister & Kiirster Il2009|) on the data. By using GLS we 



found the highest peak at P = 16.23 days with a FAP GLS = 5.5 x 1CT 6 (Fig.®. The RV variation 
has a semi-amplitude of 120 m s -1 . 



2.3. Stellar activity 

Stellar activity is an important issue when interpreting the observed RV variations. Stellar 
rotational modulation is particularly critical when probing the planetary hypothesis. It can 
mimic the star's reflex motion from a companion and so lead to a wrong interpretation of the 
RV variation. Another type of stellar activity that can produce periodic RV variation is stellar 
pulsation. We investigate both types of stellar activity of HIP 13044 in the following subsections 
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2.3.1. Rotational modulation 

A first approximation of the stellar rotation period can be obtained from the knowledge of 
the projected rotational velocity v sin i and stellar radius i?* . The maximum rotation period can be 
calculated from P TOt / sinz = 2irR*/v sini where i is the inclination angle of the stellar rotation 
axis. With a stellar radius of 6.7 R and v sin i = 10.25 km s -1 , we obtained a maximum rotation 
period of ~33 days. However, since sin 2 remains unknown, the real rotational period cannot 
be determined by this method. Obviously, a better approach should be made to determine the 
rotation period of the star. Photometric observation is a known classical method to find stellar 
rotational modulation. Periodic photometric variation due to migrating surface inhomogeneities 
(starspots) can reveal the stellar rotation. However, although there are photometric observations 
by Hipparcos and SuperWASP of HIP 13044, the observed variations, even if they are periodic, 
still cannot be directly interpreted as the result of the stellar rotational modulation. For post 
main-sequence stars like RGB and HB stars, short-period variations of hours to few days in the 
photometry are usually caused by the stellar pulsations rather than by rotational modulation. 



Bisector 



We investigated the rotational modulation by examining the variation of the line profile asymmetry 
(bisector). In this work we use the bisector velocity span (BVS) which gives the velocity 
differen ce between th e upper and lower part of the absorption line profile (for a definition of BVS 



see e.g. 



Hatzes 



1996J). We examined the BVS of the stellar spectra and searched for evidence for 
correlation between BVS and RV variation. As shown in Fig. [5] we found only a weak correlation 
between BVS and RV (correlation coefficient c = —0.13). Since the correlation coefficient is 
sensitive to the potential outliers in the data, we have carefully investigated this factor by removing 
few "outliers" in Fig. [5] The highest correlation coefficient we obtained after removing 4 outliers 
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is c = —0.4. By bootstrapping the data 100,000 times, we derived a standard deviation of 0.16 
for c. Together with the original value (c = —0.13), this confirms the weakness of the correlation 
between BVS and RV. 

Interestingly, we found that the BVS variation shows a clear periodicity of 5.02 days. The 
FAP GLS of this period is 1.4 x 10~ 5 . Fig.[6]shows the periodogram of BVS and the phase-folded 
BVS variation. Following the idea that migrating starspots on the stellar surface cause the 
variations in the line profile asymmetry or BVS, which is an effect of the rotational modulation, 
the period of Pbvs =5.02 days is a possible candidate for the stellar rotation period. In particular, 
we found no period around 16 days in the BVS variation. Hence, for HIP 13044 the period of the 
RV variation can be distinguished clearly from the anticipated stellar rotation period. 

Ca II A8498 equivalent width 

The emission cores in Ca II K (A393.4) and H (A396.7) as well as the Ca II IR triplet lines 
(A849.8, 854.2, 866.2) are well known stellar activity indicators. However, for HIP 13044 we 
found no periodic variation of the Ca II K lines. We did not use the Ca II H line for our analysis 
since it is blended by the H e line. 

Furthermore, we measured the equivalent width (EW) of Ca II IR triplets. While we did 
not find any periodicities in Ca II A854.2 and A866.2, we found periodic EW variations in Ca II 
A849.8. The LS periodogram shows a peak at 5.78 days with a marginal FAP LS of ~4% (Fig.[7J), 
whereas the GLS periodogram shows a period of 6.31 days (FAPgls = 4.4 x 10~ 6 ), as shown in 
Fig- (HI The small discprepancy between the two periods could result from the different fitting and 
weighting methods used in both periodograms. We conclude that the activity indicator Ca II 849.8 
has a period of 6.05 ± 0.3 days. 

The results from BVS and Ca II A8498 EW variation lead to the conclusion that the stellar 
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rotation period lies between 5 and 6 days. Finally, we adopted P = 5.53 ± 0.73 days, i.e., the 
mean value of the periods of both spectroscopic activity indicators as the rotation period of 
HIP 13044. 



2.3.2. Stellar pulsations 



Pulsations of Horizontal Branch stars 



Radial and non-radial stellar pulsations have been observed in post main-sequence stars, such 
as RGB, AGB stars and white dwarfs. Indeed, many HB stars are also known to be pulsators. 
However, until now there are not many studies or reported detections of RHB s tellar pulsations. 
Theoretical models of pulsations of red variable HB stars by 



Xiong et al. 



(119981) predicted that the 



oscillation periods of high-order overtones are of the order of a few tenths of a day (few hours). 



Furthermore, the oscil 
( 



Kj el dsen & Bedding 



ation frequency z/ max can be predicted from an empirical relation by 



1995). For HIP 13044 we calculated u n 



53/xHz which corresponds to 



a pulsation period of ~5.2 hours. This hypothes is has some support from the studies of HB stars 



in metal-poor globular clusters, e.g., NGC 6397 (|Stello & Gilliland 



2009). 



RR Lyrae variables 



Since the location of HIP 13044 in the Horizontal Branch is close to the RR Lyrae instability 
regime, it is also important to investiga te the pulsation s of tho se stars and compare them to the 



detected RV periodicity in HIP 13044. lHartman et al. 
in the metal-poor cluster M3. They found 



(2009) observed photometric variability 



observations, but for Ml 5 by 



Corwin et al. 



80 var iable stars with periods within 1 day. Similar 



(2008) showed that the RR Lyrae variables in Ml 5 



have pulsation periods between 0.06 and 1.44 days. These results show that a period of ~16 days 
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is not likely for oscillation characteristics of RR Lyrae stars. 



W Virginis variables 



W Virginis variable stars belong to Population II Cepheids that have pulsation periods of 10 to 30 
days. Thus, there might be a concern about the ~16 day periodicity in HIP 13044 being due to 
pulsations rather than a planetary companion. However, this requires that HIP 13044 itself must 
belong to W Virginis variables. 

Since for W Virginis stars, similar to the classical Cepheids there exist Period-Luminosity 



relation s (e.g 



period. 



Cox 



197 9b, one can calculate the expected stellar luminosity from the oscillation 



Wallerstein & Cox 



19841) listed some basic data of the Population II Cepheids in globular 



clusters. For Population II Cepheid stars with periodicities of ~16 days in the globular clusters 
M3, M10, M2 and to Cen, their absolute magnitudes My are between -2.2 and -2.7 mag. This is 
about 3 magnitudes higher than the absolute magnitude of HIP 13044 (My ~ 0.6). Less luminous 
Population II Cepheid stars with My between 0.1 and 0.6 are expected to exhibit pulsations with 
periods of 1-2 days. 



Additionaly, according to 



Pritzl et al. 



(|2003|) one can calculate the absolute magnitude of 
the star, if it indeed belongs to W Virginis variables. For P = 16.2 days, the expected absolute 
magnitude My is -1.925 mag. Because of the lack of supporting W Virginis characteristics, we 
concluded that HIP 13044 does not belong to the W Virginis variables. The 16 day periodicity is, 
therefore, not related to the stellar oscillations. 



Photometric data 



We analyzed the Hipparcos photometry data to investigate the effect of the stellar pulsations. 168 
photometric measurements of HIP 13044 are available in the public archive. The LS periodogram 
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of the Hipparcos data shows no significant signal. The highest peak corresponds to a period of 



7.1 hours with a marginal FA P=1.8%. If this perioc 



predicted pulsation frequency (|Kjeldsen & Bedding 



is tru e, it would be in agreement with the 



1995). The Hipparcos data is, however, not 



sufficient to derive more information about the stellar pulsations. 



In addition, we a 



project 



Pollacco et al. 



s o ana lyzed the photometric observations taken by the SuperWASP 



(2006). The observations were carried out in 2 blocks, each of 3-4 
months. For HIP 13044 there are 3620 photometric measurements available in the archive. After 
removing 10 outliers, probably caused by systematic errors, we used 3610 data points for the 
periodogram analysis. With such a large number of measurements, a possible 16-d periodicity in 
the photometric variations should be easily detectable. 

In Fig. [9] we present the LS -periodogram of the SuperWASP data. The LS -periodogram 
shows few marginal significant (FAP~1%) intra-day periodicities. However, it is not clear, if 
these signals can be attributed to the stellar oscillations as predicted by the theory. Interestingly, 
we observe significant peaks at 3.5 days (FAP=2 x 10~ 4 ) and 1.4 days (FAP=5 x 10~ 4 ). The two 
signals are most likely harmonic to each other. Since the FAPs of both signals are of the same 
order (10~ 4 ), the determination of the exact value of the period is still difficult. It is also possible 
that these periods correspond to oscillation overtones of very high order. Most important, however, 
is the fact that no signal at 16.2 days is detected in the LS -periodogram of the SuperWASP 
photometry data. Due to the absence of a periodic photometric variation of 16.2 days, stellar 
pulsations are not likely to be the source of the observed periodic RV variation. 



2.4. Concluding remarks 



Based on the spectroscopic and photometric analysis described above, stellar rotational 
modulation and pulsations can be ruled out as the source of the periodic RV variation. The 
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spectroscopic stellar activity indicators show a period of 5-6 days (bisectors and Ca II A849.8) 
whereas the photometric data shows periodicities of intra-day to few days, that can be clearly 
distinguished from the period of the RV variation. We conclude that the observed 16.2 days RV 
variation cannot be due to rotational modulation or stellar pulsations, but is rather caused by the 
presence of an unseen planetary companion. 



Acknowledgements 

We thank our MPIA colleagues: W. Wang, C. Brasseur, R. Lachaume, M. Zechmeister and 
D. Fedele for the spectroscopic observations with FEROS. We also thank Dr. M. Perryman, Dr. E 
Bear, Dr. N. Soker and Dr. P. Maxted for the fruitful discussion, comments and suggestions that 
helped to improve this paper. 



-22- 



Table. S2 
RV measurements of HIP 13044 



Julian Date 


RV 


a 


Julian Date 


RV 


a 


-2450000 


-300 km/s 


m/s 


-2450000 


-300 km/s 


m/s 


5103.8853 


4225.32 


50.50 


5225.6069 


4288.14 


58.04 


5106.7933 


4258.74 


47.50 


5254.5296 


4003.23 


56.92 


5164.7565 


4282.91 


50.86 


5255.5402 


4097.93 


52.17 


5168.6543 


4169.13 


50.22 


5256.5457 


4178.24 


46.06 


5170.7212 


4102.30 


50.42 


5258.5407 


4240.39 


56.84 


5172.6968 


4081.84 


63.61 


5262.5523 


4279.91 


55.03 


5174.6528 


4062.56 


42.69 


5263.5499 


4276.28 


52.83 


5175.5576 


4116.08 


51.69 


5268.5164 


4066.90 


48.92 


5197.5703 


4282.29 


47.90 


5269.5159 


4101.27 


69.23 


5199.5236 


4285.80 


48.63 


5396.8722 


4240.15 


53.09 


5199.5666 


4264.87 


44.12 


5397.8666 


4196.54 


75.11 


5200.6345 


4223.49 


51.08 


5399.8604 


4089.90 


111.8 


5202.6003 


4094.38 


46.72 


5400.9307 


4006.32 


52.99 


5203.6050 


4017.79 


54.59 


5401.8926 


4095.39 


53.14 


5204.6531 


4148.22 


64.52 


5402.8419 


4081.08 


113.26 


5223.6119 


4017.52 


54.59 


5403.8388 


4153.56 


85.38 


5224.5837 


4203.70 


51.66 


5404.8850 


4232.36 


57.91 


5225.5329 


4180.37 


50.41 


5405.8472 


4249.60 


78.77 



-23- 



0.6 



0.4 



0.2 



CD 

o 



0.6 



0.4 



0.2 



~ i — i — i — i — ii " — i 1 — i — i — i — ii 1 1 — i — i — mm 

Generalized 
Lomb-Scargle 
Periodogram of RV 




10 100 
Period (days) 



1000 



Fig. 4. — Upper panel: Generalized Lomb-Scargle (GLS) Periodogram of the RV measurements of 
HIP 13044. The window function is shown in the lower panel. The GLS periodogram shows a signifi- 
cant peak at ~16 days (dotted lines). The window function does not show any peak around this period. 
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RV-300 (km s" 1 ) 

Fig. 5. — The plot of BVS versus RV shows a weak correlation between both quantities. With all 36 
points we calculated a correlation coefficient of -0.13. A carefull rejection of 4 potential "outliers" yields a 
correlation coefficient of -0.4, which is slightly higher, but still insignificant. 
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Fig. 6. — Periodogram analysis of the BVS variation showing a periodicity of 5.02 days with a FAP of 
1.4 x 1CT 5 (upper panel). The lower panel shows a phase-folded plot of the BVS variation. The semi- 
amplitude of the variation is ~ 70ms -1 . 
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Fig. 7. — Analysis of the Call A849.8 EW variation by using a Lomb-Scargle periodogram. The LS- 
periodogram shows a peak at 5.78 days with a marginal FAP of 4%. The lower panel is the plot of the EW 
variation of Call A849.8 that is phase-folded with a period of 5.78 days. 
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Fig. 8. — Analysis of the Call A849.8 EW variation by using the Generalized Lomb-Scargle periodogram. 
The period found in the GLS -periodogram is only slightly different from that obtained from the LS- 
periodogram. The GLS -periodogram shows a peak at 6.31 days with a FAP of 4.4 x 1CT 6 . The phase-folded 
plot with P = 6.31 days is shown in the lower panel. 
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Fig. 9. — Lomb-Scargle periodogram of the SuperWASP photometry data. The periodogram shows only 
marginal significant (FAP~1%) peaks that correspond to intra-day periodicities (upper panel). It is not clear, 
if these signals are related to the stellar oscillations predicted by the theory. Nevertheless, in the lower panel 
the periodogram shows two significant peaks at 3.53 days (FAP=2 x 1CT 4 ) and 1.39 days (FAP=5 x 1CT 4 ). 
These signals are probably harmonic to each other. In this extensive data set no photometric variation with 
a period of ^16 days has been detected (dotted lines). 
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